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Abstract 

We show that if the permittivity profile of a planar dielectric medium is an analytic function in 
the upper (lower) half complex position plane then it won’t reflect radiation incident from the left 
(right), whatever the angle of incidence. Consequently, using the spatial Kramers-Kronig relations 
one can derive a real part of a permittivity profile from some given imaginary part (or vice versa), 
such that the reflection is guaranteed to be zero. This result is valid for both scalar and vector 
wave theories, and may have relevance for efficiently absorbing radiation, or reducing the reflection 
from bodies. 
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A wave propagating through an inhomogeneous medium is almost always reflected to 
some degree. This is often practically undesirable, but it is well known that in the case of an 
abrupt jump in the material parameters the reflection can be suppressed through applying 
an anti-reflection coating [TJ. However, less seems to be understood about what is required 
for a generic inhomogeneous medium not to reflect any radiation. Having said this, there are 
some famous examples of non-reflecting material profiles. One long-known example is the 
hyperbolic secant profile, which can be found in Landau and Lifshitz [2j and has been very 
clearly discussed by Lekner |3j (see [4j for an experimental realisation). More recently the 
design technique of transformation optics m has been a significant development, giving us 
a recipe for finding inhomogeneous, anisotropic materials (transformation media) that reflect 
no radiation whatever the incident field ©EE]. In the same vein perfectly matched layers [9], 
are a known family of anisotropic lossy media that are closely connected to transformation 
media, and absorb a wave without producing any reflection [TUHT3] . Another property of 
inhomogeneous media that can give an absence of reflection is PT-symmetry jlTMIfi] . This is 
a symmetry where the real and imaginary parts of the permittivity are engineered such that 
they are invariant under a simultaneous inversion of space and reversal of time. For complex 
permittivities this requires regions of gain (Im[e(a;)] < 0) as well as loss (Im[e(a:)] > 0). 
PT-symmetry guarantees zero reflection in some cases DU HU, and has been found to 
be related to the use of complex coordinates in transformation optics [T9] . Metamaterials 
allow for the realisation of such inhomogeneous permittivity and permeability profiles [ 20] 
through the use of specially designed sub-wavelength elements, and this may allow for the 
exploration of these new methods for suppressing the reflection of waves. In particular 
recent work on so-called 1 dispersion engineering ’ IZD has seen the simultaneous control of 
the real and imaginary parts of the permittivity and permeability which is necessary for the 
implementation of the following theory. 

Here we investigate the general problem of finding isotropic permittivity profiles with a 
combination of real and imaginary parts such that the reflection is zero. For planar media 
we find the very general condition that when the profile is an analytic function in the upper 
or lower half complex position plane, and therefore obeys the Kramers-Kronig relations 
in space, the reflection from respectively the left or from the right vanishes, whatever the 
angle of incidence. We note at the outset that this condition is only sufficient and not 
necessary for zero reflection. As a corollary of our finding, if the real part of such a non¬ 
reflecting permittivity profile is symmetric about some point in space then the corresponding 
imaginary part always turns out to be antisymmetric about this point, thus exhibiting PT- 
symmetry. Therefore one aspect of this work is that, similar to the findings of Castaldi et. 
al., it also points to a relationship between the use of complex coordinates, the absence of 
reflection, and PT-symmetry. 

Consider a monochromatic electromagnetic wave propagating in the x-y plane within a 
medium with an inhomogeneous permittivity e(x) that tends to a constant positive value e?, 
as x —> Too. The magnetic permeability is unity /i = 1. A schematic of this situation is 
shown in figure [l](i) . The two polarizations are TE (electric field along z) and TM (magnetic 
field along z). For the TE polarization we can write the electric field as 


E z (x,y ) = e z (x)e lkyy 

and the ^-dependent amplitude e z obeys the ID Helmholtz equation 


( 1 ) 



( 2 ) 
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FIG. 1: (i) A wave propagates in the x—y plane in a medium characterized by an inhomogeneous 
permittivity e(x) (inhonrogeneity indicated by the blue shading), (ii) Real (blue) and imaginary 
(red) parts of the permittivity profile e(x) given by equation ( [l3| ) for the parameters A = 2.0 and 
£ = 0.1A with A = 2n/ko. 


In the above equation the permittivity has the assumed form of the positive background 
contribution plus a spatially varying part 

e(x) = e b + a(x), (3) 

and the wave-number K is 

K = - kl (4) 

with ko — bj/c. The spatially varying part of the permittivity a(x) vanishes at large distances 
from the origin, where the held is made up of plane waves exp(±i/bx). 

Now suppose that we have a right-going wave that comes from infinity x = — oo and is 
incident onto the inhomogeneous permittivity profile. The effect of the permittivity profile 
is to produce a scattered held e s , and we can write the total held as 

e z (x) = E 0 e lKx + e s (x). (5) 

where K > 0. Inserting ([5]) into ([2]) we find the inhomogeneous differential equation that 
governs the scattered held 

e s (x) = -kga(x)E 0 e lKx . (6) 

One well-known way to solve equation © is to expand e s as a series 

OO 

e-fa) =XX" ) ( a 0 ( 7 ) 

71=1 

where the n th term is proportional to the n th power of a. The hrst term in this series—known 
as the Born approximation in scattering calculations—can be found immediately and is 

/ dk 

—G(k)a{k - K)e ikx , (8) 


+ K 2 + kla{x) 
ax z 


3 


















where a is the spatial Fourier transform of a{x) and G(k) is the retarded Green function 


G(k) 


1 

(K + ip) 2 - k 2 


( 9 ) 


where ij is an infinitesimal positive number. Notice that if a{k < 0) = 0 then the Born 



FIG. 2: Simulation output generated using Comsol Multiphysics. Panels (i) and (ii) show the 
absolute value of the electric field of a line source (out of the page) placed at respective positions 
x = —5A and x = +5A in the permittivity profile shown in figure [ljii). The region between the 
vertical dashed lines in (i) indicates the region plotted in figure |T](ii) . Panels (iii) and (iv) are 
for identical parameters, but we have taken only the real and imaginary parts of the permittivity 
respectively. The absence of any oscillations in panel (i) shows that the reflection is completely 
suppressed for incidence from the left, for all incident angles. 


approximation to the scattered held (|8]) is made up of only right-going waves, whatever the 
value of K (i.e. whatever the angle of incidence). This means that to first order in a(x) 
there is no backscattering from such a permittivity profile. As a first order result this is not 
all that remarkable, but through examining all the other terms in the series ([7]) we can see 
that there is actually no backscattering to any order. To prove this consider the n th term in 
the scattering series 

efdU = J^J K G(k)a(k - (A/)e“* (10) 
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This term is also made up of only right-going waves if (i) the Fourier components of the 
scattered electric field el" -1 ! are zero f or left-going waves ey n ~ l \k < 0) = 0, and (ii) the 
Fourier components of the permittivity profile are also zero for left-going waves a{k < 0) = 
0. We have already established that the first term in the series (Jr]) is made up of entirely 
right-going waves when a(k < 0) = 0, and this argument shows that every successive term 
also contains only right-going waves. There is thus zero back-scattering to all orders when 
ct(k < 0) = 0. One way to understand this result is to think that when a wave scatters 
multiple times from an object, for each scattering event there is a momentum change A k 
that occurs with an amplitude proportional to a(Ak). A permittivity profile that has only 
positive Fourier components therefore cannot convert a right-going wave to a left-going one. 

We have established that if the permittivity ([3]) is such that the Fourier transform of its 
spatial dependence is zero for k < 0, a(k < 0) = 0 then a wave incident from the left onto 
such a medium does not give rise to any reflection, whatever the angle of incidence, ft might 
appear that this argument relies on a smallness condition for a(x), but in the Supplementary 
Material we give an alternative argument that does not rely on a series expansion of the 
electric field, as well as deriving two exact solutions for propagation in such profiles that 
confirm the effect. The Supplementary Material also contains a numerical investigation 
to show that an order of magnitude increase in a(x) does not disturb the non-reflecting 
behaviour (this demonstration also shows that the real part of the permittivity can become 
negative and remain non-reflecting). 

In light of these properties, such non-reflecting permittivity profiles can be generally 
written as 

/ °° dk 

-«(fc)e ite , (11) 

which is necessarily a complex function of position. The spatial distribution of the reactive 
and dissipative parts of the material response together completely suppress reflection. To 
make use of this finding, we note that equation ( [TT| is the same in form as the relationship 
between the susceptibility in the frequency and time domains which embodies the causality 
principle [22], and one need only make the replacements k —> t and x —> oj in ( fTT| ) in order to 
recover this well known formula. As a consequence [231 123 ], the non-reflecting permittivity 
profile a(x) is an analytic function in the upper half complex position plane and satisfies 
the Kramers-Kronig relations in space 


Re[a(x)] 



7T 



s — X 


( 12 ) 


where ‘P’ indicates the principal part of the integral. Therefore if we were given some 
Im[ct(s)], say as a (square integrable) function of position, a corresponding real part can 
be constructed from (12) such that the reflection from the complex susceptibility is zero. 
We note that if the imaginary part of a(x) is symmetric about x = 0, then the real 
part calculated from (12) will be antisymmetric, and vice versa. Therefore the Kramers- 
Kronig relations generate a whole family of permittivity profiles that exhibit PT-symmetry 
(a(—x) = a*(x)). Likewise, we also have a whole family of non-reflecting profiles where 
a(x) exhibits PT- antisymmetry (a(—x) = —a*(x)), a property that has already been as¬ 
sociated with zero back-scattering in optics ESI , just as PT-symmetry [26]. Actually, even 
purely lossy periodic media can be engineered E3 such that their Bragg reflection from one 
side vanishes when the real and imaginary parts of their susceptibility are spatially out of 
phase, which is a characteristic property of Hilbert transform pairs. However, our findings 
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are more general than these known results, as they are also compatible with non-reflecting 
profiles exhibiting no definite PT-symmetry at all. 

As an initial illustration of this finding, we consider the simplest example: a permittivity 
profile with a single pole in the lower half position plane 


e(x) = e b + A 


i ~x/£ 

1 + 0/£) 2 


(13) 


where £ sets the spatial scale of the profile, and A the amplitude. Equation (13) is plotted 


in figure [T| and takes a form which would be very familiar if the x-axis represented frequency 
rather than space. The non-reflecting behaviour is demonstrated in figure [2] which shows the 
absolute value of the electric field for a point source (a line source in 2D) placed either side 


of x = 0, and compares the behaviour of the full profile (13), versus its real and imaginary 
parts separately. 

We note that, as is well known, the Helmholtz equation (J2]) is equivalent to a Schrodinger 
equation in which —a(x) plays the role of a potential profile. Thus, the spatial Kramers- 
Kronig relations generate a large family of complex non-reflecting potential profiles. Needless 
to say, as the relation given by Eq. GD is a sufficient, but not necessary condition, real non¬ 
reflecting profiles also exist which are perfectly transparent (for example, as mentioned in the 
introduction the potential V(x) = U 0 sech 2 (x/a) is known to be non-reflecting for quantum 
particles when Uq takes specific values 0Ej). 

The above analysis was carried out for TE polarization, but from our argument in terms 
of multiple scattering one might expect that these profiles are also non-reflecting for TM 
polarized waves. We now show that this is the case, given certain additional conditions on 
the permittivity. The TM polarization obeys the equation 


V ■\e- 1 (x)VH z ] + kZH z = 0 


0 


(14) 


In analogy with the foregoing argument, we write e 1 (x) = e b 1 + /3(x) and H z = h z (x)e lkyV . 
Equation (14) then takes the form 


d 2 h s 


x 


dx 2 


+ K 2 h z (x) = -e b 




(15) 


Comparing the above equation with that governing TE polarized waves ([2]) it can be shown 
that the equivalent of (10) is given by 


h ( " , (x) = € i f^f ^G(k)P(k-k')(kl + kk')h^- 1 Hk’)e‘ kx 


(16) 


Therefore if (3(k < 0) = 0 then reflection of the TM polarization is also suppressed. However, 
it is not necessarily the case that both a{k < 0) = 0 and j3{k < 0) = 0. For both equations 
to hold simultaneously we need both e(x) — e b and e~ 1 (x) — e^ 1 to be analytic functions in the 
upper half complex position plane. In particular, if e(x) satisfies the spatial Kramers-Kronig 
relations it will be free of zeros in the upper half plane when lm[e(a;)] takes only one sign along 
the real axis (a proof of this property of analytic functions can be found in [28]). Therefore 
lossy media obeying the Kramers-Kronig relations in space will not reflect radiation of 
either polarization for any angle of incidence, which establishes the generality of the finding 
illustrated in figure [3] Meanwhile for profiles exhibiting a combination of loss and gain where 
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there are zeros in the upper half position plane will suppress reflection for only one of the 
two polarizations. Note that because the TM polarization is sensitive to zeros in e(x) then, 
unlike the TE polarization, the non-reflecting behaviour is sensitive to to the value of e^. 

As an example that demonstrates the generality of our finding we now take a permittivity 
profile with an imaginary part given by 


Im[e(x)] = 


h(L — x ) 


[1 + erf(x/{)] [1 + erf ((L - x)/()] 


(17) 


which represents a smoothed triangle function, where ‘erf’ is an error function, h is the height, 
L the length, and £ characterizes the smoothness of the corners. Numerically calculating 


the integral (12) we obtain the real part of the permittivity that, when added to i times (17) 
is necessary to reduce the reflection to zero. The full function is shown in figure [3^ v), and 
unlike (13) has no definite parity symmetry. Figures [3](i-iv) then show that the resulting 


profile is non-reflecting for TM polarized waves incident from the left, but that there is 
reflection from the right and for the real and imaginary parts of a(x) taken separately. 

A general result in the theory of reflection [29] is that when waves are incident onto a 
generic planar interface, at angles close to grazing (k y ~ k 0 ) the reflectivity usually ap¬ 
proaches unity. The above findings at first sight appear to contradict this result. There 
are some other somewhat surprising features of these prohles which are related to this. For 
instance if we could construct a medium with a permittivity profile that obeys the spatial 
Kramers-Kronig relations over all frequencies (this property does not obviously contradict 
the Kramers-Kronig relations in frequency), then it would be non-reflecting for all angles 
of incidence and all frequencies. 

This surprising behaviour stems from the fact that strictly speaking the prohles we have 
calculated are of infinite extent so that there is no ‘interface’ to speak of, and no natural 
length scale associated with the prohles. The importance of the infinite extent of the prohles 
is indicated by the presence of the long tails evident in the real parts of e{x) shown in figures[I] 
and [3j In practice we must confine these infinite prohles to a finite region of space through 
truncating these tails, which can be achieved through multiplying the prohle by an envelope 
function U(x) 


a'(x) = a(x)U(x ) 
r°° (jy 

a'(k) = / —U{k' - fc)a(fc') 

Jo 27r 


where U(x) = 0 at some distance from the centre of the prohle. This truncation nat¬ 
urally introduces a length scale into the permittivity prohle, and through doing this we 
find that a'{k < 0) 7 ^ 0. However, if U(Ak) is sharply peaked around A k = 0, a' 
will only be non-zero for negative k of a small magnitude and nearly all of the non- 
rehecting behaviour can be retained. For example if we use U(x) = exp(— (x/a) 2 ), then 
U(Ak) = (a/ 2 v / 7 r) exp(—a 2 (A/c) 2 /4), which rapidly goes to zero beyond around Ak = —2/a. 
The consequence of this is that waves close to grazing will now be rehected by the prohle. 
The Supplementary Material contains a further discussion of this effect, where it is numeri¬ 


cally demonstrated that it is possible to confine the prohle (13) to a slab of a few wavelengths 


thickness, introducing reflection at close to grazing incidence, while otherwise retaining the 
non-rehecting behaviour. 

While the non-rehecting property of specihc classes of real potential prohles have 
long been studied [30] . more recently the analogous behaviour of PT-symmetric complex 
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FIG. 3: Panel (i) shows the absolute value of the field of a dipole source aligned along the y axis, 
and placed at x = —5A in front of the permittivity profile shown in panel (iii). Panel (ii) the 
same situation, but with the dipole placed at x = 5A. Panels (iii) and (iv) show the behaviour 
when we retain only the real and imaginary parts of a(x) respectively. The imaginary part of the 
permittivity profile is given by © with = 2, L = 1.5A, £ = 0.1A, and h = 1.0. The real part 
of the profile was numerically calculated using (12). The scale on the axes of the field plots is the 
same as in figure [2] 


susceptibility profiles have been considered in optics. Yet, to the best of our knowledge the 
simple and general relation here discussed between the one-sided absence of reflection and 



















the analytic extension of the spatially dependent susceptibility to one half of the complex 
position plane has not been pointed out before. We have shown how the corresponding 
Kramers-Kronig relations in space can be used to generate a large family of non-reflecting 
profiles as they provide a sufficient condition for being non-reflecting on one side. If 
the profile is also free of zeros in the upper or lower half complex position plane (the 
half plane being determined by whether reflection vanishes from the left or the right) 
then the profile is also non-reflecting for both polarizations. In practice the catch is 
that the profiles have very long tails which must be truncated, and where one chooses 
to perform the truncation determines the range of angles and frequencies that are not 
reflected. Nevertheless, the advantage of this method is that it requires us only to be 
able to manipulate the real and imaginary parts of an isotropic permittivity, and is in 
principle valid for any wave equation, including the Schrodinger equation. While the 
Kramers-Kronig relations in the frequency domain are a cornerstone of optics, it is hoped 
that the spatial Kramers-Kronig relations will provide at least some guidance and in¬ 
sight in the development of metamaterials based on judiciously chosen susceptibility profiles. 
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